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Abstract 

The 4>^ distribution of the Z/y* — >■ £'^£~ production in hadron collisions is simulated using 
a leading-order event generator, GR@PPA. The initial-state parton shower, which simulates the 
multiple QCD-radiation efltects in the initial state, plays the dominant role in this simulation. The 
simulation in the default setting agrees with the high-statistics measurement by ATLAS at LHC 
with the precision at the level of 5%. The observed systematic deviation, which can be attributed 
to the effects of ignored higher-order contributions, can be reduced by adjusting the arbitrary 
energy scales in the simulation. The agreement at the level of 1% can be achieved over a very 
wide range without introducing any modification in the implemented naive leading-logarithmic 
parton shower. 



1 Introduction 

Z boson production in hadron collisions is a unique place for studying the initial-state strong interac- 
tion (QCD) activity, because the production kinematics can in principle be unambiguously measured 
if the production is tagged by a lepton pair (electrons or muons) from the decay. The spectrum of the 
Z-boson transverse momentum, pt{Z), with respect to the beam direction, which can be reconstructed 
from the measured decay-lepton momenta, is the most important property for probing the underlying 
QCD interactions. We expect that the pt{Z) spectrum at high pt{Z) (> mz) can be well described 
by the production of a high-p-r jet (light quark or gluon) in association with the Z production, where 
mz represents the invariant mass of the Z boson. However, as pxiZ) decreases, the contribution of 
the multiple QCD radiation becomes significant and alters the naive prediction from the Z + 1 jet 
production. Finally, at pt{Z) < 10 GeV/c, the divergent Z + 1 ^et cross section is suppressed to 
form a peak structure at several GeV/c. A certain non-perturbative contribution is considered to be 
important around the peak. These features must basically be common to all interactions in hadron 
collisions since the multiple-radiation effects can be factorized. 

Monte Carlo (MC) event generators are indispensable tools for evaluating the detection efhciency 
and acceptance of experiments at high-energy hadron colliders. Since the events in the low-p^ peak 
region dominate the total yield, it is important for MC event generators to precisely reproduce the 
spectrum around the peak. MC event generators simulate the multiple radiation effects with parton 
showers (PS) based on the leading-logarithmic (LL) approximation, and non-perturbative effects with 
certain models. Thus, their performance is not trivial, particularly around the peak. The Z-boson 
Pt spectrum is a good test bench for the verification of implemented approximations and models. 
However, because of the presence of a finite error in the momentum measurement of high-momentum 
decay leptons, it is difficult for experiments to make precise measurements of the pt{Z) spectrum at 
low pt{Z). 

Alternatively, an angular variable has been introduced for probing the transverse activities in 
Z-boson production pj. The introduced variable </>* is defined as 

(/.; =tan(?iacop/2)sinr. (1) 



1 



The acoplanarity angle (/)acop is the supplement of the opening angle between the two decay leptons 
in the projection onto the plane perpendicular to the beam direction. The angle 9* is defined by the 
pseudorapidity, r/ = — ln[tan(0/2)], of the leptons (^^) as cos6'* = tanh[(r7+ — r/^)/2]. Although, in 
general, the momentum resolution deteriorates as the pt of charged particles increases, the angular 
resolution usually remains constant irrespective of pr- A transverse recoil of Z bosons more or less 
produces a finite acoplanarity between the decay leptons, and the amount of the acoplanarity depends 
on the decay angle in the Z-boson rest frame. The variable (/)* has been proposed to maximize the 
sensitivity to small Z-boson recoils, and it has an approximate correlation topT(-^) as 4>* « pT{Z)/mz- 

Subsequent to the proposal of this angular variable, the DO experiment at FNAL Tevatron, which 
provides proton-antiproton collisions at a center-of-mass energy {\/s) of 1.96 TeV, published their 
measurement on the 0* spectrum 2 . They found that the resummation calculation by RESBOS 
provides a good prediction with accuracy better than 10% over almost the entire measurement range. 
However, they also found a substantial systematic deviation between their measurement and RESBOS. 
Although such observations are interesting to be cross-checked, DO has presented the results based 
on a lepton momentum measurement that is not fully corrected for the final-state photon radiation 
(QED FSR) effects in the Z decay. Hence, an appropriate correction is required in the simulations to 
be compared. This fact causes a difficulty in the simulation side. 

Recently, a new result on the 0* spectrum with significantly higher statistics has been published 
by the ATLAS experiment [1]. Their measurement is based on proton- proton collision data at -y/s — 
7 TeV provided by CERN LHC. In contrast to DO, ATLAS has presented the results for several 
definitions of the lepton momenta. Among them, the Born-level definition, in which all QED FSR 
effects are corrected for in the data analysis, allows a straightforward comparison with simulations. 
This definition also allows to combine the results for the Z e+e^ and Z fJ-^ tJ-~ channels. Using 
the combined high-precision result, ATLAS has confirmed the systematic deviation from RESBOS 
that has been observed by DO with more stringent significance. 

In this article, we study the performance of the GR@PPA event generator by comparing its pre- 
diction on the 0* spectrum with the ATLAS measurement. GR@PPA [5J[3] is an event generator for 
multi-body production processes in hadron collisions based on the GRACE system [7]. The version 2.8 
package (U [S] supports various weak-boson production processes in which associate one-jet production 
is combined by using a jet-matching method (TUl El [H] to reproduce the weak-boson kinematics in 
the entire phase space. The package includes custom-made PS programs in order to ensure the per- 
formance of the jet matching. We have demonstrated in previous studies [HI [13] that the GR@PPA 
predictions reproduce the pt{Z) spectra measured at Tevatron and LHC with good precision without 
any tuning. 

According to the motivation to introduce the 0* variable, we focus on the spectrum in non-hard 
regions, pt{Z) < mz, roughly corresponding to (/)* < 1.0. The spectrum is predominantly determined 
by the initial-state PS simulation in this region. The multiple QCD-radiation effects simulated by PS 
are evaluated by resummation calculations [3J [TJ] on the basis of next-to-leading logarithmic (NLL) 
or further higher order approximations. Although the theoretical basis of practical PS programs 
is still limited to the LL order, higher-order contributions are partially imported in terms of the 
angular ordering in some PS programs [151 116j . Hence, a question arises: is it necessary to include 
the effects beyond the LL order also in PS simulations for hadron collision interactions? The answer 
depends on the accuracy that we require. Before answering this question, we have to investigate how 
accurately the LL approximation can reproduce actual phenomena. Since the primitive features of 
the LL approximation are strictly preserved in the GR@PPA PS, the present study will provide us 
with an opportunity to probe the capability of the LL approximation itself. 

The GR@PPA PS covers the energy region at > (5 GeV)^. Perturbative and non-perturbative 
contributions at smaller are simulated down to the hadron level by feeding the generated event 
information to the general-purpose event generator PYTHIA [TB]. We have found in a previous 
study [12] that the additional PS simulation in PYTHIA is not effective in the pt{Z) spectrum, 
whereas the simulation of non-perturbative effects is effective around the pt{Z) peak; it alters the 
peak position. Hence, the verification of the PYTHIA model for non-perturbative effects, in the 
combination with the GR@PPA PS, is also an important subject in the present study. 

The rest of this article is organized as follows. The simulation using the GR@PPA event generator 
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is described in Sec. [2l The result in the default setting is compared with the ATLAS measurement. 
Possible alterations in the simulation to achieve better matching with the measurement are discussed 
in Sec. [31 and the discussions are concluded in Sec. HI 

2 GR@PPA simulation 

In the present study, we use the 2.8.4 update of the GR@PPA event generator, in which the im- 
provements established in the study of diphoton production |17| are migrated to the previous 2.8.3 
update [S]. However, the impact of this migration is negligible for the Z-boson production kinemat- 
ics, because the improvements are mainly concerned with the final-state PS. The event generation is 
carried out for the 7-TeV LHC condition, proton-proton collisions at -y/s — 7 TeV. The Z + jet and 
Z + 1 jet production processes are combined by using the jet-matching method. The decay of the Z 
boson is included in the matrix elements for the event generation, and all Z bosons are assumed to 
decay to an electron pair. The photon exchange contribution and its interference with the Z exchange 
are also included in the matrix elements. The events are generated without any cut, except for the 
constraint on the invariant mass of Z bosons, 66 < mz < 116 GeV/c^. 

MRST20071omod [ig in the LHAPDF 5.8.4 ^ library is used for the parton distribution func- 
tion (PDF) in protons, as in the previous study [13]. The choice of PDF is not important in the 
simulation of the transverse QCD activity, except for the spectrum at very high pT (-^ fnz) where 
the Z + 1 jet contribution becomes dominant. The initial-state QCD radiations are simulated by 
using the backward-evolution PS (QCDPSb) included in the GR@PPA package. The energy scales, 
the factorization scale {fJ^p) and the renormalization scale (/ifl), are chosen as iip = A*-R = '^^z- The 
energy scale of the initial-state PS is always equal to fj,F in our jet- matching method. The final-state 
PS (QCDPSf) is also activated. Although it is not important in the present study, the final-state PS 
energy scale is also set to mz- 

The generated events are exported to PYTHIA 6.425 (TB] in order to add simulations at small 
along with simulations of hadronization and decays. The default setting in PYTHIA is unchanged 
except for the setting of PARP(67) = 1.0 and PARP(71) = 1.0, as in the previous studies [TB]. 
The ATLAS measurement is presented for the quantities of Born-level leptons, in which the QED 
FSR effect in the Z decay is corrected for. In order to follow this definition, we look for an electron 
pair from a Z-boson decay before the QED FSR simulation in the PYTHIA event record. The pair 
is the same one as that generated by GR@PPA, but their rest frame is boosted and rotated by the 
simulations in PYTHIA. Since these electrons are defined at the Born level, they can be replaced with 
muons. Hence, we refer to these electrons as generic leptons {£^) in the following discussions. 

The event selection is applied using the momenta of these leptons according to the definition of 
ATLAS; that is, pr > 20 GeV/c and jr/j < 2.4 are required for both leptons. The invariant mass is 
already constrained in the event generation. The selected events are binned according to the binning 
on (p* defined by ATLAS to derive the {1 / a)da / del)* spectrum. 

Figure [1] shows the ratio of the ATLAS measurement with respect to the GR@PPA simulation 
result. The results for the Z — >■ e"'"e~ and Z — >■ pJ^ fJi~ channels are combined in the ATLAS measure- 
ment. We can see the agreement fairly better than 10% between the simulation and measurement over 
the entire measurement range. However, we also see an apparent systematic deviation between them. 
The tendency is similar to that observed in the comparison with RESBOS in the ATLAS paper [i], 
although the deviation is slightly larger in the present result. 

The most remarkable feature of the result in Fig. [1] is that the measurement/simulation ratio is 
significantly smaller than unity at medium 0*, 0.1 ^ (j)^"^ 0.5. This means that the event fraction in 
the medium pt{Z) range, 10 < pt{Z) ^ 50 GeV/c, is too large in the simulation. The fraction be- 
comes smaller than the measurement in the peak region, 0* < 0.1, to compensate for the enhancement 
since the compared distributions are normalized. We found in a previous study [12[ that the pt{Z) 
spectrum at pt{Z) > 10 GeV/c is not altered by the PYTHIA simulation. Hence, the observation 
implies that the PS radiation effect is too strong in the GR@PPA event generation. 

Figure [5] shows the result of the comparison in three rapidity regions of the lepton pair. The 
rapidity is defined by the total energy E and the momentum along the beam direction pz as y = 
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Figure 1: {l/a)da/d(j)* spectrum normalized to the prediction from the GRQPPA simulation with the 
default setting. The result for the spectrum measured by ATLAS is plotted. The solid and dotted 
bars illustrate the predictions from GR@PPA with extreme settings of the arbitrary energy scales: 
solid bars for = ^fi = 0.5mz, and dotted bars for = ^ji = 2.Qmz- 
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(a) \y\ <0.8 

(b) 0.8< |y| <1.6 

(c) |y| >1.6 



Figure 2: {1 / a)da / dc/)* spectrum normalized to the prediction from the GR@PPA simulation with the 

default setting. The result for the ATLAS measurement is separately plotted in three rapidity regions 
of the lepton pair: (a) < 0.8, (b) 0.8 < < 1.6, and (c) \yu\ > 1.6. 
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(1/2) ln[{E + pz) / {E — pz)]. We can see that the overall behavior is almost independent of the rapidity, 
except in the large 0* region, (j)* > 1.0, where the Z + 1 jet contribution is dominant. The gluon 
density in PDF contributes to the Z + 1 jet cross section together with quark densities, whereas the 
total yield is predominantly determined by quark densities. Thus, the fraction in the large cj)* region 
should be altered by the change in the gluon/quark ratio in PDF in the relevant parton momentum 
range. The observation at large (p* would provide additional information on PDF. However, we do 
not discuss further details of this issue in this article. Detailed discussions should be carried out using 
the pt{,Z) distribution. 

It is to be noted that the error bars illustrated in Figs. [T] and [5] show the quadratic sum of all errors: 
the statistical and systematic errors of the measurement and the statistical error of the simulation. 
The simulation is based on 10^ events generated by GRQPPA, and the number of events is reduced to 
about 5 X 10^ after the event selection. Because the ATLAS measurement is based on a huge amount 
of data, approximately 3 x 10^ lepton-pair events, the contribution of the simulation statistics is not 
negligible in the illustrated errors. However, this contribution does not affect the above discussions 
since the observed systematic deviation is considerably larger than the estimated error. 

3 Optimization of the simulation 

The Z production cross section evaluated on the basis of the lowest-order qq Z interaction neces- 
sarily has a dependence on the factorization scale {(J-f), since PDFs are provided as a function oi fip. 
In principle, this scale can be chosen arbitrarily. This dependence is mostly canceled by the inclusion 
of Z + 1 jet production in our jet-matching method [10115]. and the inclusion of the Z -f 1 jet extends 
the pt{Z) spectrum to high pt{Z) to cover the entire phase space. We can achieve good matching 
between the Z -I- jet and Z -I- 1 jet in the pt{Z) spectrum by introducing an appropriate parton 
branch model [11] [12]. However, because the cancellation is realized only at the first order of the 
QCD coupling, as, a substantial dependence remains not only in the total cross section but also in 
the Pt{Z) spectrum [THHl]. 

The evaluated cross section also has an ambiguity due to the arbitrariness in the renormalization 
scale (/ii?). The as value in matrix-element (ME) calculations depends on ^r. Since — >■ Z is an 
electroweak interaction, the ambiguity affects the Z production cross section only through the Z + 1 
jet production. Thus, the effect is small but substantial. 

These ambiguities depending on the choice of the energy scales, jip and /i/?, are caused by the 
lack of higher-order contributions. The variation in the simulation result reflects the uncertainty in 
the theory that the simulation is based on. In our jet- matching method, we take /i/j = /i^ in order 
to equalize the radiation probability in the PS applied to the Z -I- jet and in the MEs for the Z + 
1 jet at their boundary, = fi^. The energy scales are usually set to the "typical" energy of the 
interaction of interest, because non-collinear contributions ignored in PS and PDF become significant 
and higher-order contributions are minimized around this energy. Hence, we take fin ^ ~ "niz as 
the default for Z production, as described in the previous section. 

Theoretical uncertainties are usually evaluated by multiplying the default choice of the energy 
scales by factors of 0.5 and 2.0. The results for these extreme choices, ^ir = hf = O.Sm^ and 
^J■R = t^F = 2.0mz, are illustrated in Fig. [T]using solid and dashed bars, respectively. We can see that 
the ATLAS measurement lies within the range of the uncertainty represented by these two results in 
the range (j)* <1.0. 

The theoretical uncertainty is rather large because the MEs are evaluated only at the tree level. 
However, the main purpose of MC event generators is not to provide theoretical predictions, but to 
provide a tool for analyses in experiments. The simulations are desired to reproduce actual measure- 
ments as precisely as possible in this usage. Adjustments of parameters in the simulation should be 
allowed to satisfy this requirement. Of course, the adjustments have to be carried out in a reasonable 
range or within the range that the underlying theory allows. Looking at the results in Fig. [TJ we 
speculate that better matching may be achieved with energy scales slightly smaller than the default 
values, although it is not trivial that the observed systematic deviation vanishes with such a simple 
adjustment. 
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Figure 3: {1 / a)da / dff)^ spectrum normalized to the prediction from the GR@PPA simulation with 
the adjusted energy scales, = /i^ = O.TSm^. Two results for the ATLAS measurement are 
plotted. Filled circles show the ratio to the simulation employing the default primordial fc^ simulation, 
kr = 2.0 GeV, in PYTHIA, while open circles show the ratio to the simulation with fcy = 2.4 GeV. 
Dotted lines indicate the ±1% deviation between the measurement and simulation. 
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filled circle: kj = 2.0 GeV (default) 
open circle: /f,- = 2.4 GeV 

{a)|y|<0.8 

(b) 0.8 < |y| <1.6 

(c) |y| >1.6 



Figure 4: {l/a)da/d(j)* spectrum normalized to the prediction from the GR@PPA simulation with the 
adjusted energy scales, ^ir = ^j,f = 0-75mz- The notation is the same as that in Fig. [31 but the results 
are separately shown in three rapidity regions of the lepton pair: (a) < 0.8, (b) 0.8 < \yu\ < 1.6, 
and (c) \yff \ > 1.6. Dotted lines indicate the ±1% deviation between the measurement and simulation. 
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We have tried several settings and found that we can obtain an almost optimum result with the 
choice of ^ir = hf = 0.75mz- The result with this setting is plotted in Fig. [3] with filled circles. In 
this result, the simulation statistics are improved by a factor of five with respect to the study in the 
previous section. The resultant statistical error is approximately 1/3 of that from the measurement; 
thus, it is almost negligible. We can see that the dip that is observed in Fig. [T] in the medium-0* 
region, 0.1 <(/>*< 0.5, has vanished. The simulation matches the measurement with accuracy better 
than 2% through the measurement range up to (p* « 1.0. Note that the vertical scale is expanded in 
this plot with respect to Fig.[T] Figure S] shows the results in the three rapidity (yu) ranges. We can 
see that the adjustment has reasonably improved the matching in all yu bins. 

In spite of the observed overall improvement, the measurement/simulation ratio still shows a 
substantial systematic discrepancy from unity in the small-t/)* region, (/)* < 0.1. This small-0* region 
corresponds to the peak region in the priZ) spectrum, pt{Z) < 10 GeV, where simulations for non- 
perturbative effects in PYTHIA based on certain models are effective. Among various simulations 
applied in PYTHIA, the primordial-fcr simulation plays the most important role for determining the 
Pt{Z) spectrum in this region |12) . This simulation is considered to include non-perturbative effects 
together with small-Q^ perturbative effects below the PYTHIA-PS cutoff, = (1.0 GeV)^. The 
applied model is very simple; an additional pr is assumed for each parton in protons according to the 
Gaussian distribution having the standard deviation of fcy- The default setting is kr — 2.0 GeV. 

The result in Fig. [3] shown with filled circles indicates that the PYTHIA model for the primordial 
kx is almost appropriate. However, there is still some room for optimizing the kx parameter. The 
open circles in Fig.[3]show the result when we increase the fcy value to 2.4 GeV. The plots are slightly 
shifted to the left in order to ensure that the changes are visible. Although the matching deteriorates 
at very small </>* (< 0.01), where even the angular measurement may be ambiguous, a substantial 
improvement can be seen in the region of 0.01 <</>*< 0.1. As a result, the measurement/simulation 
ratio distribution becomes flat within ±1% over a very wide range, 0.01 < (/>^ ^ 1.0. The results in 
three separate yu regions after the kx adjustment are also illustrated with open circles in Fig. |4l We 
can see substantial improvements that are similar to the improvement observed in Fig. [3] 

4 Conclusion 

The measurement of the angular variable 4>*, which is defined for the production of Z bosons in 
hadron collisions decaying to a lepton pair, provides us with rich information concerning the initial- 
state QCD activities, especially at small pt{Z) where the pt{Z) measurement is ambiguous. It 
has been demonstrated in this article that the normalized 4>* spectrum, {1/ ij)da /d(j)^, measured by 
the ATLAS experiment at LHC can be reproduced by the simulation based on the GR@PPA event 
generator with the precision better than 5% over almost the entire measurement range where the 
spectrum is predominantly determined by multiple QCD radiations simulated by PS. 

Simulations based on MC event generators are desired to reproduce measurement data as precisely 
as possible when they are used as tools for measurements. We have shown that the observed disagree- 
ment between the measurement and simulation is well within the theoretical uncertainty, which is 
represented by the variation due to the change of arbitrary energy scales in the event generation, 
and we have demonstrated that the agreement can be improved to the level of 1% by adjusting 
these scales. It should be emphasized that this level of precision has been achieved with a primitive 
leading-logarithmic PS implemented in the GR@PPA event generator. 

We use the PYTHIA event generator for simulating non-perturbative and perturbative effects 
at small (< (5 GeV)^) in our simulation. The successful simulation at small </)* implies that 
the simulation in PYTHIA is basically appropriate. However, the remaining fine structure in the 
measurement/simulation ratio suggests the necessity of a fine tuning of the kx parameter in the 
primordial-fcT simulation in PYTHIA. The change of the parameter kx to 2.4 GeV from its default 
value of 2.0 GeV improves the matching in the relevant 0* range. 

In this article, we have discussed the capability of the GR@PPA event generator as a tool for 
Z-boson production measurements. We anticipate that a similar level of precision can be achieved for 
other processes, although the optimum choice of the parameters would be dependent on the process 
because the higher-order contributions would be different. The collision-energy dependence of the kx 
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parameter may also be an interesting subject for future investigation. The GR@PPA event generator 
is being developed to provide a tool for such applications and studies. The most up-to-date version 
of the package is available from the Web pag^ll- 
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